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ABSTRACT 

We report the detection of WASP-35b, a planet transiting a metal-poor 
([Fe/H] = —0.15) star in the Southern hemisphere, WASP-48b, an inflated planet 
which may have spun- up its slightly evolved host star of 1.75Rq in the Northern hemi- 
sphere, and the independent discovery of HAT-P-30b / WASP-51b, a new planet in 
the Northern hemisphere. Using WASP, RISE, FTS and TRAPPIST photometry, with 
CORALIE, SOPHIE and NOT spectroscopy, we determine that WASP-35b has a mass 
of 0.72 ± 0.06M/ and radius of 1.32 ± 0.03i?j, and orbits with a period of 3.16 days, 
WASP-48b has a mass of 0.98±0.09Mj, radius of L67iO.08.Rj and orbits in 2.14 days, 
while WASP-51b, with an orbital period of 2.81 days, is found to have a mass of 
0.76 ± 0.05M/ and radius of 1.42 ± 0.04i?j, agreeing with values of 0.71 ± 0.03M/ 
and 1.34 ± 0.07Rj reported for HAT-P-30b. 
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Subject headings: Stars: planetary systems 



1. Introduction 



Around 120 transiting exoplanets have now been discovered^, having both their radii and 
masses known, a nd may be subjected to further investigations into their bulk and atmospheric com- 
posit ion (see e.g. lCharbonneau et al, booi 120051 : iFortnev et aklhoOTl : lHaswellll2010l : ISeager fc Demingl 
20101 ) . They show a diverse range of densities and thus internal compositions. Over 80% of those 
with masses greater than O.lMj have radii greater than that o f Jupiter, and many have extremely 



low d ensities which even co reless models struggle t o explain (IFortnev et al 



2007 



Burrows et al 



2011 



and 



20071 ). for e xample TrES-4b (IMandushev et alJl2007l >. WASP-17b (|Anderson et alJl201CL 
Kepler- 7b (jLatham et al.ll20ld ). More transiting exoplanets are vital to constrain the formation 
and evolution of such planets. 

Several possible effects on planetary radii have been discussed in the l iterature, inc l uding 



Guillot & Showman 


2002). 


ohmic heating from the ( 


flows (Batvgin & Stevenson 


2010: 


Batvgin et al. 


2011 


of close-in exoplanets (Bodenheimer et al. 


2001, 


2003 



1996 



inflation due to strong h eating by irradiation received from the host star (jGuillot et al 

coupling of magnetic fields and atmospheric 
and tidal heating d ue to the circularisation 



Jackson et al.ll2008h . 



Most of the known transiting exoplanets orbit very close to their host star, which produces 
strong tidal forces betwee n them. T he tidal interactions may result in orbital circularisation, 
synchronisation and decay (jPontl l2009l ) . An apparent relationship between period and eccentricity 
of non-transiting exoplanets, such that planets very close to their host stars have circular orbits 
while those farther from the stars show a wi de range of e ccentricities, seems to provide evidence for 
the tidal circularisation of planetary orbits (|Mazehll2008l ). However, it had generally been assumed 
that the masses of pl anets are too small to synchronise the stellar rotation with the planetary orbit 
through tidal forces (jMazehl 120081 ) . Some possible exceptions exist for planets with large masses, 
for example the r Boo system wit h a 4Mj planet, for which there is evidence for stellar rotation 
synchronisation (IFares et al.l 120091). R ecently. CoRoT-llb has been announced to have a 'peculiar 



tidal evolution' (jGandolfi et al.ll201ll ). where the stellar rotation period of 1.7 days is now shorter 
than the planetary orbital period of 3 days, due to the strong tidal forces causing orbital decay of 
the planet from a previously approximately synchronised orbital period. The possibility of perhaps 
unknown close-in exoplanets affecting the stellar spin rate of their host stars produces uncertainty 
in the use of the gyro chronological age estimate of a star, based on the normal spin down of a star 
as it ages. 

One motivation of the WASP project is to discover enough transiting exoplanets with a wide 
range of orbital and compositional parameters to allow analyses that may distinguish between 
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different models. Recent WA SP discoveries hel p fill in areas of planetary mass-radius paramter 



space, for example WASP-39b (jFaedi et al,ll201ll ) is the least dense sub-Saturn mass planet known. 



Here, we announce WASP-35b, WASP-48b and WASP-51b (= HAT-P-30): the positions and 
magnitudes of each star are given in Table [H We set out the photometric and spectroscopic ob- 
servations obtained for each system, before presenting the resulting stellar and system parameters. 
WASP-48 appears to be a slightly evolved star with a large radius of 1.75Rq and a very short 
rotation period of around 7 days for its age, around 8 Gyr. WASP-48b may therefore have spun-up 
its host star despite a mass of only around lMj, due to the star's inflated radius, although the 
system has not yet reached synchronisation of stellar rotation period with planetary orbital period 
(2.15 days). As the host star continues to evolve, increasing significantly in radius, it may engulf 
WASP-48b, which is currently orbiting at a distance of around 4 times the stellar radius, unless 



orbit al decay due to tidal interactions, such as may have happened to CoRoT-llb (jGandolfi et al 
boill ). push the planet far enough away to survive. 



Table 1: Star locations. 



Star 


RA (J2000) 


Dec (J2000) 


V mag (NOMAD) 


WASP-35 


05 04 19.6 


-06 13 47 


10.95 


WASP-48 


19 24 39.0 


+55 28 23 


11.66 


WASP-51 


08 15 48.0 


+05 50 12 


10.36 



Observations 



2.1. Photometric Obervations 



The WASP-North and South observatories are located at La Palma, Canary Islands and SAAO, 
South Africa respectively, and each consist of eight 11cm telescopes of 7.8° x 7.8° field of view each, 
on a single fork mount. The cam eras scan repea t edly t hrough eight to ten sets of fields each night, 



taking 30 second exposures. See iPollacco et al.l (120061) for more detai ls on the WAS P 



pro 



ect and 

the data reduction procedure, and ICollier Cameron et al.l (|2007l ) and IPollacco et al.l (|2008l ) for an 
explanation of the candidate selection process. 

High precision transit light curves of WASP-35, WASP-48 a nd WASP-51 were obtained with 



the R ISE instrument mounted in the 2.0m Liverpool Telescope ([Steele et al.ll2008l ; I Gibson et al 



20081 ) at the Observatorio del Roque de los M uchachos on La Pal ma, Canary Islands. The data 
we re redu c ed wi th the ULTRACAM pipeline (jDhillon et al.ll2007l ) following the same procedure 



MS 



Barrosl ( 2011 



Further high quality photometry of WASP-35 was obtained with the 0.6m 
TRAnsiting Planets and Planeteslmals Small Telescope (TRAPPIST) at La Silla Observatory, 
Chile, and with the LCOGT 2.0m Faulkes Telescope South (FTS) at Siding Spring Observatory. 
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Photometric observations for each of WASP-35, WASP-48 and WASP-51 are detailed in Table 
[2l and shown in Figures [I] to [3j 

Table 2: Photometric Observations 



Instrument 


Date(s) 


Num datapoints 


WASP-35 


WASP-N (2 cams) 
WASP-S (1 cam) 


10/2008-01/2009 
10/2009-01/2010 
10/2009-01/2010 


16,748 


RISE (V+R) 


30/11/2010 


630 


TRAPPIST (I+z) 


23/12/2010 


742 


FTS (z) 


01/01/2011 


205 


TRAPPIST (I+z) 


11/01/2011 


300 


WASP-48 


WASP-N (3 cams) 


05/2007-07/2007 
05/2008-08/2008 
05/2009-06/2009 
05/2010-09/2010 


56,036 


RISE (V+R) 


01/07/2010 


770 


WASP-51 


WASP-N (1 cam) 


11/2009-03/2010 


4,453 


RISE (V+R) 


12/01/2011 


1,245 



2.2. Spectroscopic Observations 

We obtained 3 radial velocity measurements of WASP-35 on the 25 and 26 January 2010 
with the Fibre-fed Echelle Spectrograph (FIES) on the 2.5m Nordic Optical Telescop^]. FIES 
was used in medium resolution mode (R=46,000) with simultaneous ThAr calibration, and the 
observations were reduced using the FIEStool packag^f], then cross-correlated with a high signal-to- 
noise spectrum of the Sun to obtain the radial velocities. A further nine datapoints were obtained 



2 operated on the island of La Palma jointly by Denmark, Finland, Iceland, Norway and Sweden, in the Spanish 
Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias 

3 http://www.not. iac.es/mstruments/fies/fiestool/FlEStool. html 
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Fig. 1. — Upper plot: WASP-35 discovery lightcurve folded on the orbital period of P = 3.16 d. 
Middle plot: WASP data binned in phase, with bin-width = 11 minutes. Lower plot (top to 
bottom): transit obtained with the RISE instrument on 30 November 2010, TRAPPIST on 23 
December 2010 and 11 January 2011, and FTS on 1 January 2011, lightcurves offset for clarity. 
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Fig. 2. — Upper plot: WASP-48 discovery lightcurve folded on the orbital period of P = 2.14 d. 
Middle plot: WASP data binned in phase, with bin-width = 8 minutes. Lower plot: transit obtained 
with the RISE instrument on 1 July 2010. 
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Fig. 3. — Upper plot: WASP-51 discovery lightcurve folded on the orbital period of P = 2.81 d. 
Middle plot: WASP data binned in phase, with bin-width = 10 minutes. Lower plot: transit 
obtained with the RISE instrument on 12 January 2011. 
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with the CORALIE echelle spectrograph mounted on the Swiss telescope at the ESO La Silla 
Observatory in Chile on 5 January 2010 and between 6 and 14 February 20 10 and the resulting 
spectra were proc essed with the standard data reduction C ORALIE pipe line (jBaranne et al.lll99a : 
Mayor et al.ll2009l ). plus a correction for the blaze function (|Triaudll2010l ). 



WASP-48 was observed spectroscopically with the SOPHIE fibre-fed echelle spectrograph 
mounted on the 1.9 3m telescope at Observatoire de Haute-Provence (OHP) (jPerruchot et alj|2008l : 
Bouchv et al.ll2009l ). Fourteen datapoints were taken in High Efficienty mode (R = 40,000) between 
19 April and 18 October 2010, obtaining a signal-to-noise ratio of around 35 for each. Wavelength 
calibration with a Thorium- Argon lamp was performed to allow the interpolation of the spectral 
drift of SOPHIE. 

Radial velocity observations of WASP-51 were made with the CORALIE and SOPHIE spec- 
trographs. Nine datapoints were taken with SOPHIE between 16 October and 25 November 2010, 
using the high resolution mode (R = 75,000) and obtaining a signal-to-noise ratio of around 20 
for each measurement. Ten further datapoints were obtained with CORALIE between 7 December 
2010 and 4 January 2011. 

The CORALIE guiding camera revealed a faint star located at about 1.5" from WASP-51 which 
was observed simultaneously with WASP-51 through the 3" fiber of the SOPHIE spectrograph. 
Thus, we performed radial velocity blend simulations (Santerne et al, in preparation) in order 
to investigate whether small bisector variations seen in SOPHIE data are due to this secondary 
blending star, though the bisector spans are compatible with no variation within the 2a level. This 
simulation consisted of producing SOPHIE cross-correlation functions (CCFs) for WASP-51 with 
the unblended orbital solution seen by CORALIE and blending a single solar star with different 
systemic radial velocity, distance from WASP-51 and vsini. With 10,000 Monte Carlo simulations 
in the ranges 34.6-54.6 kms -1 for the systemic RV, 1-5 times the distance of WASP-51 and 4- 
40 kms -1 for the vsini of the blending star, we didn't find any significant model that reproduces 
the bisector effect seen in the SOPHIE data. For the case of WASP-51, the effect of the blending 
star on the bisector is less than a few m.s -1 . We concluded that the faint star 1.5" away didn't 
affect the SOPHIE CCF of WASP-51. 

The radial velocities measured using these spectra are given in Tables [3l 0] and [5l and shown 
as a function of orbital phase in Figures HJ [5] and El along with radial velocity residuals and bisector 
span variations. The radial velocity measurements show low-amplitude variations, compatible with 
the existence of orbiting plane ts, and are not co rrelated with the bisector spans (available for 
CORALIE and SOPHIE data) (IQueloz et al.ll200ll ). which rules out stellar blend configurations. 
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Table 3: Radial velocity measurements of WASP-35 



BJD-2 400 000 


RV 


^ rv,V 


BS 


Instrument 




(km s — ) 


(km s — ) 


(km s — ) 




5201 7298 


17.8174 


0085 


-0 0225 


CORALIE 


5221.5256 


17.7386 


0.0096 




FIES 


5222.3276 


17.6211 


0.0079 




FIES 


5222.5131 


17.6318 


0.0093 




FIES 


5233.6480 


17.8050 


0.0133 


-0.0020 


CORALIE 


5235.6597 


17.6897 


0.0084 


-0.0446 


CORALIE 


5236.5322 


17.8342 


0.0078 


-0.0531 


CORALIE 


5237.5430 


17.7133 


0.0081 


-0.0172 


CORALIE 


5238.5912 


17.6439 


0.0080 


0.0069 


CORALIE 


5239.6245 


17.8183 


0.0088 


-0.0402 


CORALIE 


5240.5331 


17.7754 


0.0128 


-0.0087 


CORALIE 


5241.5318 


17.6446 


0.0075 


-0.0201 


CORALIE 



Table 4: Radial velocity measurements of WASP-48 taken with the SOPHIE spectrograph 



BJD-2 400 000 RV ctrv BS 







(km 


s" 1 ) 


(km s 1 ) 


(km s 1 


5305 


,6441 


-19 


,711 


0.024 


-0.029 


5323 


,5842 


-19 


,648 


0.021 


0.015 


5328 


,5784 


-19 


831 


0.025 


-0.083 


5409 


4890 


-19 


,667 


0.026 


-0.071 


5427, 


,5503 


-19 


,747 


0.023 


-0.062 


5429, 


,4215 


-19 


,809 


0.026 


-0.044 


5431, 


,4254 


-19 


831 


0.033 


-0.014 


5432 


3967 


-19 


,589 


0.036 


-0.059 


5476, 


,2714 


-19 


,745 


0.021 


-0.002 


5477, 


,2679 


-19 


591 


0.026 


0.044 


5482, 


2945 


-19 


,594 


0.022 


-0.029 


5484, 


,2801 


-19 


,529 


0.020 


0.042 


5485, 


,2880 


-19 


,823 


0.026 


-0.027 


5488, 


,2831 


-19 


,569 


0.036 


0.009 
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Table 5: Radial velocity measurements of WASP-51 



BJD-2 400 000 


RV 


0"RV 


BS 


Instrument 






(km s _1 ) 


(km s _1 ) 


(km s _1 ) 




5485. 


,6249 


44.611 


0.018 


-0.023 


SOPHIE 


5489, 


6119 


44.783 


0.018 


-0.017 


SOPHIE 


5495 


6696 


44.708 


0.015 


0.100 


SOPHIE 


5513. 


,6757 


44.596 


0.013 


-0.013 


SOPHIE 


5519. 


,6870 


44.657 


0.014 


-0.002 


SOPHIE 


5522 


6691 


44.677 


0.015 


0.038 


SOPHIE 


5523 


,7176 


44.754 


0.023 


0.084 


SOPHIE 


5524, 


,6337 


44.593 


0.022 


-0.055 


SOPHIE 


5525 


,6785 


44.771 


0.021 


-0.003 


SOPHIE 


5527 


6506 


44.586 


0.015 


0.048 


SOPHIE 


5537 


,8023 


44.7111 


0.0062 


-0.0046 


CORALIE 


5538 


,8006 


44.5639 


0.0061 


-0.0118 


CORALIE 


5542 


,8449 


44.7469 


0.0073 


0.0240 


CORALIE 


5544 


7669 


44.6278 


0.0065 


0.0087 


CORALIE 


5545 


,8487 


44.7761 


0.0066 


0.0016 


CORALIE 


5547, 


,8576 


44.6832 


0.0074 


-0.0027 


CORALIE 


5561 


,8052 


44.6526 


0.0072 


0.0241 


CORALIE 


5562 


,8426 


44.7309 


0.0070 


-0.0010 


CORALIE 


5563 


,7958 


44.5945 


0.0083 


0.0155 


CORALIE 


5565, 


,7984 


44.7353 


0.0070 


0.0170 


CORALIE 
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Fig. 4. — Radial velocity measurements for WASP-35: red, open diamonds denote NOT obser- 
vations, blue circles are from CORALIE. The solid line shown on the top plot is the best-fitting 
MCMC solution with eccentricity fixed to zero. The centre-of-mass velocity, 7 = 17.71 km s , was 
subtracted. The middle plot shows the bisector span variations against radial velocity measure- 
ments. The uncertainty in the bisector span measurements are taken to be twice the uncertainty in 
the radial velocity measurements. The lower plot shows the radial velocity residuals against time. 
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Fig. 5. — Radial velocity measurements for WASP-48: details are as for WASP-35, with centre-of- 
mass velocity, 7 = -19.68 km s" 1 here. 
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Fig. 6. — Radial velocity measurements for WASP-51: red diamonds and blue circles denote OHP 
and CORALIE observations, respectively. Details are as for WASP-35, with centrre-of-mass veloc- 
ity, 7 = 44.68 km s _1 here. 
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Results and Discussion 



3.1. Stellar Parameters 



The individual NOT/FIES spectra of WASP-35 were co-added to produce a single spectrum 
with an average S/N of around 100:1. For WASP-48, the SOPHIE spectra were co-added to produce 
a spectrum with typical S/N of around 80:1. For WASP-51, the CORALIE spectra were co-added 
to produce a single spectrum with an average S/N of around 150:1. The standard pipeline reduction 
products were used in the analysis. 



The analysis was performed using the methods given in lGillon et al.l (|2009i ) . The Ha line was 
used to determine the effective temperature, T e g, while the Na I D and Mg I b lines were used 
as surface gravity, logg diagnostics. The elemental abundances were determined from equivalent 
wid th measuremen t s of se veral clean and unblended lines. Atomic line data was mainly taken from 
the iKurucz fe Belli (119951) comp ilation, but with updat ed van der Waals broad e ning coefficients 



for lines in 



Barklem et al 



|2000j) and log gf values from iGonzalez fe Laws! (J2000|) , iGonzalez et al 



( 2001 ) or Santos et al.l (2004 ). A value for microturbulence, £tj was determined from Fe I using the 



method of iMagainl ()1984l ) . The parameters obtained from the analyses are listed in Table [UJ The 
quoted error estimates include that given by the uncertainties in T e g, logg and £tj as well as the 
scatter due to measurement and atomic data uncertainties. 

The projected stellar rotation velocity (v sini) was determined by fitting the profiles of several 
unblended Fe I lines. Values for macroturbulence, w ma cj of 3.3 ± 0.3 kms -1 , 3.2 ± 0.3 kins" 1 and 
4.0 ± 0.3 kms " 1 were assumed fo r WASP-35, WASP-48 and WASP-51, respectively, based on the 
calibration bv lBruntt etall (hold ), and instrumental FWHM of 0.13 ± 0.01 A, 0.15 ± 0.01 A and 
0.11 dh 0.01 A, respectively, determined from the telluric lines around 6300 A. Best fitting values 
of vsini= 3.9 ± 0.4 kms -1 , 12.2 ± 0.7 kms -1 and 3.6 ± 0.4 kms -1 , respectively, were obtained. 



WASP-51b has been recently announced as HAT-P-30b (jjohnson et al.l l201ll ). A discrep- 
ancy between their results an d ours is the yalue of host star metallicity: we find a value of 
[Fe/H] = -0.08 ±0.08 whereas Ijohnson et al.l (hoilh repo rt fFe/Hl = 0.13 ±0.08, from an analysis 



using SME (Spectroscopy Made Easy. IValenti fe Piskunovl ()1996l )). The difference is partly due to a 
differ ent solar metall i city y alue adopted in the spectral analysis, where w e use log A(Fe) = 7.54±0.03 



which is 0.04 dex higher than the current lAsp 



from iBiemont et al.l (|1991 

value of log A(Fe) = 7.50 ± 0.04. We continue to use the IBiemont et al 



und et, al.1 (120091 ) solar 



(]199ll ) value for consis- 



tency between analyses rather than adjusting to new values as they are announced. Additionally, 
we allow the microturbulence value to float and be fitted, findin g a best-fit value of 1.3 kms -1 , in 
excellent agreement with t he calibration in iBruntt et al.l ([20101 ), compared to t he as sumed value 
of 0.85 kms -1 adopted by IValenti fe Fischeii (j2005l ) and used in I Johnson et al.l ([20 111 ). Using the 
lower value of microturbulence would increase our [F e/H] abundance by 0.07. The slight differences 
in T e g and log g values, where I Johnson et al.l (120111 ) find T e g = 6300 ± 90 and log g = 4.36 ± 0.03, 
result in a difference of a further 0.03 dex. These three factors combine to reduce the discrepancy 



to only 0.07 dex, which is within the quoted uncertainties. 

Table 6: Stellar parameters of WASP-35, WASP-48 and WASP-51 from Spectroscopic Analysis. 



Parameter 



WASP-35 



WASP-48 



WASP-51 



T e S 

log g 
v s'mi 



6050 ± 100 K 
4.4 ± 0.1 

1.3 ± 0.1 kms^ 1 

2.4 ± 0.6 kms" 1 



6000 ± 150 K 
4.5 ± 0.15 
1.3 ± 0.2 kms- 1 
12.2 ± 0.7 kms" 1 



6250 ± 100 K 
4.3 ± 0.1 
1.3 ± 0.1 kms" 1 
3.6 ± 0.4 kms" 1 



[Fe/H] 


-0.15 


± 


0.09 


-0.12 ± 0.12 


-0.08 ± 0.08 


[Na/H] 


-0.16 


± 


0.11 


-0.14 ± 0.07 


-0.11 ± 0.05 


[Mg/H] 


-0.01 


± 


0.11 


-0.09 ± 0.07 


-0.08 ± 0.08 


[Al/H] 








0.02 ± 0.14 




[Si/H] 


-0.06 


± 


0.05 


0.04 ± 0.08 


-0.01 ± 0.09 


[Ca/H] 


-0.01 


± 


0.12 


-0.10 ± 0.11 


0.05 ± 0.12 


[Sc/H] 


-0.02 


± 


0.07 


0.13 ± 0.09 


0.01 ± 0.06 


[Ti/H] 


-0.03 


± 


0.10 


0.02 ± 0.13 


-0.03 ± 0.08 


[V/H] 


-0.13 


± 


0.06 


-0.17 ± 0.16 


-0.10 ± 0.09 


[Cr/H] 


-0.10 


± 


0.09 


0.09 ± 0.16 


-0.01 ± 0.06 


[Mn/H] 


-0.29 


± 


0.08 




-0.27 ± 0.11 


[Co/H] 


-0.22 


± 


0.10 






[Ni/H] 


-0.16 


± 


0.07 


-0.18 ± 0.09 


-0.10 ± 0.06 


log A(Li) 


2.35 ± 0.08 


< 1.3 


2.96 ± 0.08 


Mass 


1.10 ± 0.08 


1.09 ± 0.08 M 


1.20 ± 0.09 M 


Radius 


1.09 ± 0.14 


1.09 ± 0.14 R Q 


1.28 ± 0.17 R Q 



Note: Mass and radius estimates using the iTorres et al.1 (120101 ) calibration 



The lithium abunda nce in WASP-35 implies an age of >2 Gyr according to the calibration of 

23 ±6 . 5 days 



Sestito fc Randichl (|2005l ) while the measured v sin % gives a rotational period of P, 



rot 



assum ing i = 90°, which yields gyrochronological age of 



fi +5 



Gyr using the relation of 



Barnes 



(|2007l ) . A lack of stellar activity is indicated by the absence of Ca II H+K emission in the spectra. 



The non-detection of lithium in the WASP-48 spectrum suggests that the star is several Gyr 
old, and the lack of any Ca H+K emission is consistent with this. However, the stellar rotation 
rate of 4.5 ± .6 days from the -usini measurement, implies an age of only ~ 0.2^'^ Gyr using the 
Barnes! (|2007l ) gyrochronological relation. This discrepency is discussed further below. 



The presence of strong lithium absorption in the WASP-51 spectrum suggests that the star 
is <1 Gyr old. The stellar rotation period of 18.0 ± 3.1 days found from the vsini measurement 
results in an age of ~ 5.6jVg Gyr, though this is can be regarded as an upper limit, since the stellar 
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rotation rate may be higher due to the unknown value of stellar inclination (i). 



3.2. System Parameters 

We simultaneously analysed all photom etry and radial velocity data using a Markov-Chain 



Monte Carlo (MCMC) analysis, as set out in ICollier Cameron et al.1 (|2007l ). modified to calculate 



the m ass of the host star with a calibration on log p and [Fe/H] as described in lEnoch et al 



(|2010l ). The photometry provides the stellar density while the effective temperature and metallicity 



of the star are obtained through spectral analysis, given above. 

We initially allowed the value of the eccentricity of each planetary orbit to float, which resulted 
in values of e = 0.057+°;°^, e = 0.058±g;g|| and e = 0.040+°;°^ for WASP-35b, WASP-48b and 
WASP-51b, respectively . The eccentricity of all three orbits are consistent with zero based on 



the Lucy-Sweeney test (ILucy h Sweeney! Il97ll ) on the data for each system, which yields 100% 
probability that these eccentricity values could be due to chance if the underlying orbit were circular 
for all three systems. With no strong evidence for eccentricity in the orbits, we then performed 
second analyses of all systems, fixing the eccentricities to zero. 

The MCMC results for WASP-48b reveal a slightly evolved star with a large radius of 1.75 Rq 
(discussed below). This radius is around O.QRq greater than expected for that of a star of the same 
mass on the Main Sequence. To explore this, we performed a further analysis, imposing a Main 
Sequence constraint on WASP-48, such that a Gaussian prior is used on the stellar radius at each 
step in the analysis, centred on radius R* = M^' 8 , again allowing the eccentricity to float. This 
forced the stellar radius estimate down to 1.13i?0, artificially inflating the eccentricity estimate to 
0.38 ± 0.04 to compensate, since the stellar radius is calculated from the photometry, including a 
term 1 ^ S1 ^ 1 : an eccentricity of 0.4 produces a stellar radius value 0.65 times that of the estimate 
with eccentricity at zero. However, the Lucy-Sweeney test result is still 100%, and there is no 
evidence for such a high eccentricity in the radial velocity curve, where the x 2 values are 33.3 and 
8.6 for eccentric and circular orbits, respectively. Imposing the Main Sequence constraint with 
the eccentricity held fixed at zero decreases the stellar radius to 1.31 ± O.O9i?0 and decreases the 
impact parameter from 0.7 to 0.4 to compensate: such a reduction in impact parameter leads 
directly to a stellar radius estimate 0.70 times smaller. However, the much lower impact parameter 
estimate produces a model that does not well match the high quality follow up photometric data: 
compare the fit to ingress and egress in Figure El showing the lower impact parameter model, with 
that in Figure [21 Since neither analysis with Main Sequence constraint imposed produces realistic 
models for the photometry and spectroscopy, the results from the analysis with no Main Sequence 
constraint are clearly preferred. 

We also tested the imposition of the Main Sequence constraint on WASP-35 and WASP- 
51, producing negligible change in all parameters, well within the parameter uncertainties, which 
implies that the photometry is sufficient to well constrain the results. 
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Fig. 7. — RISE photometry shown with the low impact paramter (b=0.4) model; a poor fit com- 
pared to the fit with b=0.7 in Figure [21 



The best-fit parameters for each systems, with eccentricity fixed to zero and no Main Sequence 
constraint imposed, are given in Table WASP-35 is found to have a mass of 1.07 ± 0.02 M 
and radius of 1.09 ± 0.02 R Q , while WASP-35b has a mass of 0.72 ± 0.06 Mj and radius of 
1.32 ± 0.03 Rj, giving a density of 0.32 ± 0.03 pj. For WASP-48, the best-fit model gives a 
stellar mass of 1.19±O.O4M and radius of 1.75±0.07H Q , giving a stellar density of 0.22 ± 0.02p Q . 
WASP-48b is found to have a mass of 0.98 ± 0.09Mj and radius of 1.67 ± 0.08Rj, giving a density 
of 0.2li&,. 

The best-fit result for WASP-51 gave a stellar mass 1.18±O.O3M , stellar radius 1.33±O.O3i? , 
planetary mass 0.76 ± 0.05Mj, planetary radius 1.42 ± 0.04i?j, giving a density of 0.26 ± 0.02pj. 
Under the name HAT-P-30b, based on HAT-5 (Arizona) and HAT-9 (Hawaii) photometry of a total 
of around 3,200 good datapoints, follow-up high quality photometric observations wi th Kepler-Cam 



on the FLWO 1.2m telescope, and Keck and Subaru radial velocity measurements, iJohnson et al 



(|201ll ) report a planetary mass of 0.71 ± 0.03Mj, and radius of 1.34 ± 0.07Rj, which agree with 
our values. 



3.3. Isochrone Analysis 



The values for stellar density, effective temperature and metallicity from the MCMC analyses 
(MS constraint off, eccentricity fixed to zero) were used in an interpolation of the Padova stellar 



Table 7: System Parameters from simultaneous MCMC analysis of all photometric and spectroscopic data. 



1—* O T* O YY1 AT QT 


oy iiiuui 






vv -o i 


Period (days) 


T) 

r 


o.lbitwb ± U.UUUUUz 


Z.14obo4 ± U.UUUUUa 


Z.olUbUo ± U.UUUUUo 


Iransit r^pocn (HJDj 


T 


5531.47907 ± 0.00015 


53b4. 55043 ± 0.00027 


5571.70057 ± 0.00020 


Iransit duration (days) 


D 


n -i 0^70 i n nnn c 

0.1278 ± 0.0005 


0.1327 ± 0.0013 


n nnon i n nnno 

0.0920 ± 0.0008 


Planet/ Star area ratio 




n m r - a i n nnni 

U.U154 ± U.UUUl 


n nnn/ 1 i n nnno 

U.UU9o ± U.UUUz 


n m oo i n nnno 

(J.(Jlz2 ± (J.UUUz 


Impact Parameter 





U.oU ± U.(J4 


U. (6 ± U.Uo 


U.o / ± U.U1 


btellar Keflex Velocity (ms ) 




94 - 82 -6.92 


13b.02I 11 44 


96.70_ 604 


Centre-of-mass Velocity (kms" 1 ) 


7 


l/./z ± U.U1 


i a coo _i n nm 
— 19. boo ± U.UU1 


44.0/bo ± U.UUUl 


Orbital separation (AU) 


a 


0.04317 ± 0.00028 


n no aaa i n nnn /i o 

U.Uo444 ± 0.1)1)04,3 


n n /f i i o i n nnno 1 

0.04118 ± 0.00031 


Urbital inclination (deg) 


i 


87.96I 025 


on nn+0.69 
80 - 09 -0.55 


do /IO+0.14 

82.48I 017 


Urbital eccentricity 


e 


U (fixed) 


(fixed) 


U (fixed) 


Orbital distance / radius ratio 


a/it* 


o c o t n i n 

8.5o ± 0.19 


yi oq+0.24 
4 - 23 -0.19 


C7 1 n 1 *7 

b.b7 ± 0.17 


Stellar mass (Mq) 


M* 


1 r\7 _i_ n no 
1.U7 ± 0.02 


1 m i n r\A 

1.19 ± U.U4 


1.18 ± 0.03 


Stellar radius (Rq) 


■ft* 


i no _i_ n no 
i.uy zt u.uz 


i 7^ -i- n n.7 

1.IO x U.U / 


i -i- n n^ 


Stellar surface gravity (log g & ) 


log 5* 


4.40 ±0.01 


4.03 ± 0.03 


4.26 ±0.01 


Stellar density (p & ) 


P* 


0.83 ± 0.03 


0.22 ±0.02 


0.50 ±0.02 


Stellar metallicity 


[Fe/H] 


-0.15 ±0.09 


-0.12 ±0.12 


-0.08 ±0.08 


Stellar effective temperature 


T eff 


5990 ± 80 


5920 ± 150 


6250 ± 100 


Planet mass (Mj) 


M p 


0.72 ±0.06 


0.98 ± 0.09 


0.76 ±0.05 


Planet radius (Rj) 


R p 


1.32 ±0.03 


1.67 ±0.08 


1.42 ±0.04 


Planet surface gravity (log gj) 


log % 


2.98 ± 0.04 


4.03 ± 0.03 


2.93 ±0.03 


Planet density (pj) 


Pp 


0.32 ± 0.03 


0.21 ±0.04 


0.26 ±0.02 


Planet temperature (A=0, F=l) (K) 


Teq 


1450 ± 20 


2030 ± 70 


1710 ± 30 
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evolution tracks (jGirardi et al.1 12002| ; iMarigo et al.l 120081 ) , shown in Figures [U [9] and [TUJ This 
resulted in an age for WASP-35 of 5.01 ± 1.16 Gyr and a mass of 1.03 ± 0.04M Q . 

WASP-48 was found to have an age of 7.9±f;g Gyr, and mass of 1.08 ± 0.06Mq, supporting the 
results from the MCMC analyses and lack of lithium and Ca H+K that WASP-48 is an old, slightly 
evolved star. However, the inferred age from the stellar rotation period is still only 0.6^q' 2 Gyr when 
based on a rotational period of 7.2 ± 0.5 days found from the larger stellar radius of 1.75 ± 0.07-Rq 
(with v sini= 12.2 ±0.7 kms -1 , as before). This apparent contradiction may be due to stellar spin- 
up by the planet. The close orbits of most known transiting exoplanets produce strong tidal forces 
between the planet and host star, which lead eventually to orbital circularisation, synchronisation 
and decay. There is evidence for tidal circularisation among the known transiting exoplanets, 
seen in the period-eccentricity relation where planets with very short orbital periods have circular 
orbits while those with longer periods show a range of eccentricities, but the synchronisation of 
stellar rotation is generally ruled out for most known planet-star systems, with timescales orders of 



magnitude larger than the Hubble time (IMazeh 



2008 



Pontll2009l ) . A few exceptions exist where the 



companion is massive, e.g. the unique case of r Boo for which the planetary orbital period coincides 
closely with the mean rotation period of the ste llar differential rotation pattern determined from 
the Zeeman-Doppler imaging (jFares et al.l l2009l ). WASP-48b orbits with a period of 2.14 days 
(a = 0.034 AU), and while it is not a massive companion, the evolved nature of the star giving a 
radius of 1.75 Rq results in a relatively short stellar tidal s ynchron i sation ti me of lljj Gyr. A full 



tidal analysis will be discussed in the forthcoming paper of iBrownl (jin prepl ). 



For WASP-51, the interpolation yields an age of 4.01J-2 Gyr, and mass of 1.16 ± 0.04Af ( 



z . 



An alternative analysis was also performed using the Yonsei-Yale isochrones (jDemarque et al 



20041 ). which produced an age estimate of 4^ Gyr and a mass estimate of 1.04 ± 0.04Af Q for 



WASP-35, 6.0±?;g Gyr and 1.14 
51. 



+0.06 
-0.07 



M for WASP-48 and 3.0lJ;g Gyr and 1.20t^M Q for WASP- 



+0.03 ; 



The results of both isochrone analyses are consistent within uncertainties, though the Padova 
interpolations yield higher ages for each system than the Yonsei-Yale interpolations, and mass 
values agree with the calibrated values from the MCMC analysis of 1.08 ± O.O3M , 1.19 ± 0.04M Q 
and 1.18 ± O.O2M for WASP-35, WASP-48 and WASP-51, respectively. 



4. Summary 

We have presented observations of WASP-35b, WASP-48b and WASP-51b, new transiting 
exoplanets of 0.72 Mj and 1.32 Rj, 0.98 Mj and 1.67 Rj, and 0.76 Mj and 1.42 Rj, respectively. 
We have obtained photometry from both WASP-North and -South instruments, as well as high- 
quality follow-up photometry from RISE, FTS and TRAPPIST. Radial velocity measurements from 
the FIES, SOPHIE and CORALIE instruments allowed us to confirm the planetary nature of the 
companions. 
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Fig. 8. — Isochrone tracks from lMarigo et al.l (j2008l ) for WASP-35. Isochrones (solid lines, left to 
right) are: 1.0, 1.99, 3.16, 5.01, 6.30, 7.94, 10.0 Gyr. Evolutionary tracks (dashed lines) are for 1.2, 
1.1, 1.0, 0.9 M Q stars. 



WASP-48b may have spun-up its slightly evolved host star, which shows a high v sin i value, 
otherwise contradicting the estimated stellar age of around 8 Gyr. WASP-51b has been recently 
announced as HAT-P-30b, with mass 0.71M/ and radius 1.34i?j, values that agree with our results, 
within uncertainties. 

WASP-North is hosted by the Isaac Newton Group on La Palma and WASP-South is hosted by the South 
African Astronomical Observatory (SAAO) and we are grateful for their ongoing support and assistance. 
Funding for WASP comes from consortium universities and from the UK's Science and Technology Facilities 
Council. The RISE instrument mounted on the Liverpool Telescope was designed and built with resources 
made available from Queen's University Belfast, Liverpool John Moores University and the University of 
Manchester. The Liverpool Telescope is operated on the island of La Palma by Liverpool John Moores 
University in the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias 
with financial support from the UK Science and Technology Facilities Council. We thank Tom Marsh for the 
use of the ULTRACAM pipeline. TRAPPIST is a project funded by the Belgian Fund for Scientific Research 
(FNRS) under the grant FRFC 2. 5. 594. 09. F, with the participation of the Swiss National Science Fundation 
(SNF). M. Gillon and E. Jehin are FNRS Research Associates. The research leading to these results has 
received funding from the European Community's Seventh Framework Programme (FP7/2007-2013) under 
grant agreement number RG226604 (OPTICON). 
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Fig. 9. — Isochrone tracks from iMarigo et al.l ( 20081 ) for the slightly evolved WASP-48 (~ 8 Gyr, 
1.75^0). Isochrones (solid lines, left to right) are 3.16, 3.98, 5.01, 6.30, 7.94, 10.0, 12.5 Gyr. 
Evolutionary tracks (dashed lines) are for 1.2, 1.1, 1.0, 0.9 Mq stars. 
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